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Abstract: This paper studies the critical exponents of a doubly degenerate parabolic equation with 
nonlinear boundary flux. The equation can be used to describe the nonstationary flow 
in a porous medium of fluids with a power dependence of the tangential stress on the 
velocity of the displacement under polytropic conditions. This study is thus meaningful 
in various branches of applied mathematics. We obtained the global existence exponent 
po and critical Fujita exponent pe by constructing various self-similar supersolutions and 
subsolutions of the objective doubly degenerate parabolie equation. The main results of 
this paper are as follows the equation exists global solutions if 0 < p < po; the equation does 
not exist nontrivial and nonnegative global solutions if po < p < pc; the equation exists 
global solutions for small initial valus while does not exists global solutions otherwise if 
P > Pe- 
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1 Introduction 
In the paper, we consider the critical exponents to the doubly degenerate parabolic equation 
Ut = (juslP(u™ Ja), (x,t) € Rt x (0,T), (1) 


with 
—luz|?(uw™). (0, t) = uP(0,t), tel, T), 


u(x, 0) = uo(z), xeERt, 


where the parameters m > 1, p, 8 > 0. 

The equation (1) can be used to describe the nonstationary flow in a porous medium of 
fluids with a power dependence of the tangential stress on the velocity of the displacement 
under polytropic conditions, appearing in several branches of applied mathematics. 

It is well known that the critical exponent problem was first studied by Fujital!] in 1966 for 
Cauchy problem of the semilinear equation up = Au + u”. The critical exponent of the doubly 
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degenerate equation 
ut = div(|Vul?Vu™) +u”, (x,t) € S= R” x (0,00), (3) 


with 8 > 0, m > 1, p> 0 was proved in [2,3] that pe = m + 8 + B+? In [4] it was determined 
that the degenerate problem 


Ut = (lus us), (z, t) eRt x (0,7), (4) 


with nonlinear boundary —|uz|?u,(0,t) = u?(0,t), t € [0, T) and initial data u(x, 0) = uo(z), x 
€e R+, where 8, p > 0 admits the critical exponents po = 2042 and pe = 28 + 2, which can be 
obtained by taking m = 1 in (1). A similar work was arad in a [5]. The critical exponents for the 
heat equation with a mixed nonlinear Dirichlet-Neumann boundary condition were considered 
in [6]. We refer to [7] for a survey on studies of critical exponents in particular!®1 for the 
results of critical exponents to degenerate parabolic equations or systems. 

In this paper we will prove that the critical exponents of (1) should be 


2B6+m+1 
Po = B}? 
Together with such po and pe, the main results of this paper are stated by the following theorems. 
Theorem 1 Any solutions of (1) are global if 0 < p < po. 
Theorem 2 Any nontrivial solutions of (1) blow up in finite time if po < p < De. 
Theorem 3 The solutions of (1) are global with small initial data and non-global with 
large initial data if p > pe. 


Pe = 28-+m +1. 


2 Proof of the main results 


We prove the main results of the paper in this section. 
Proof of Theorem 1 Construct 


L 
T(x, t) = e (M +e hee") = eltem, 


with 
M = max {|luoll™, 1}, L= (M +1) FED man, 
seq ee pete is (Btm), (5) 
A simple computation shows 
T > ke*t M”, (mf (@™) x), < (B+ yee ? Pk- Dt+(mk— 2ljt Pl -1 (6) 
with m > 1, and hence 
T > (\te|9(U™)z),, in Rt x (0,7) (7) 


by (5). Moreover, with p < po and (5), it holds that 


—|%z|P(a™).(0,t) > GP(0,t), for te (0,7) (8) 
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on the boundary, and 
U(x,0) > uo(xz), on Rt (9) 


for the initial data. We conclude from (7)-(9) that T is a global supersolution of (1) whenever 
O<p<po. 
On the contrary, we have the following lemma. 
Lemma 1 The solutions of (1) blow up in finite time with large initial data if p > po. 
Proof Construct u(x,t) = (T — t) ESE), € = 2(T - t), where f is the compactly 
supported function to be determined with 
= B+1 be p-B-m 
(B+ 2)p — (28 +m-+1)’ (8 +2)p— (28 +m+1) 


It follows from (10) that 








k (10) 


k+1=A(k+l) +km42l, B(k+1l)+km+l=pk. 
After a direct computation, we have 
uy < (ltte|8(u"™)a)g» (2t) ERt x OT), —[uelP(u™)o(O,t) < uP(0,t), t€0,7) (11) 
if f(€) satisfies 
KFO HEFE SUF PUY EO, -EET O < FO): 


Set f(€) = A(cA9+m~?) — €)4 = Az} with 


B+1 _ mab+? 


= Bama S ko 





Q 


and A > 0 to be determined later. Noticing (a — 1)+ ma — 1 = a and z < cAPtm-—1, we have 
KFE) + lef’ E) — (IF PUF™))' (E) < 27 AP (ke — mah) < 0. 


If we choose l 
A> (mab +! i P)e) FEF=CRFRED , 


then —(|f’|9(f™)')(0) — fP (0) < 0. Thus, u(x,t) is a non-global subsolution of (1) with the 
initial data large that uo(x) > u(x,0) on Rt. 

In the sequel, we follow the techniques in [4] to prove Theorems 2 and 3. 

Proof of Theorem 2 Construct 


up(z,t)=(tt7)°9(O, gE) =C- E) E=a(Q +7) (12) 


with 7, c > 0 and 


1 B+1 B+2 À —6-1 Fai 
= 2 = L = ——,, = = C= —(L i 3 
A m+28+1’ B+m-1 y B+i (ah 7) ) 


It is easy to verify that g(€) satisfies 


rg(E) + Ag (E) ~ (I'P YY E = 9, —(I9’1? gy) ©) = 9? (0). (13) 
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Since g’(0) = 0, the self-similar solution ug(z,t) satisfies ug(0,t) = 0. Without loss of gener- 
ality, suppose u(0, to) > 0 for some to > 0. Therefore, there exist large 7 > 0 and small c > 0 
such that u(z,to) > ug(z,to) for all z > 0. A direct computation shows that up(z,t) is a 
subsolution of (1) on R+ x (to, +00), namely, u(x,t) > ua(z,t), (x,t) € Rt x (to, +00). 

We claim that there exist t* > to and T large enough such that ugp(z,t*) > u(x,0) on Rt. 
In fact, wp(z,t*) > u(x,0) is true provided that 


(t* 47) FET > TO OTR GED, (th +T) ET K T BRED (14) 





Since p < pe = 28 + m + 1 implies 


B+1 X p-B-—m 
(B+1)p—(26+m+1)~ (6@4+1)p—(26+m+1)’ 


there exist t* > to and T large enough that (14) are both satisfied. We conclude that u(x, t*) > 
up(a,t") > u(z,0), x € Rt. Based on Theorem 1, it can be conducted that every nontrivial 
nonnegative solution of (1) blows up in finite time if po < p < De. 

Proof of Theorem 3 Construct G(x, t) = (t+7)7*h(€), € = x(t +7)~!, where constants 
T > 0, and k, l are defined as (10). If h(€) satisfies 


-kh — leh! > (IWEK), (WIP (kY (0) > hP(0), (15) 
then @ is a supersolution of (1) with small initial uo(x). We set 
h(€) = A(C((da)? — (€ + a)7){) = Ag(€ + a) (16) 
with C, g, y, L defined as (12). From (13), we know that g(€ + a) satisfies 


(I'P) = —Ag— AlE + a)g', —(|9/P(g™)') (a) = g(a), 


2 (17) 
g'(E + a) = —CLy(é + a)! (da)? — (€ + a)7) $7. 
By a simple computation, we know that the first inequality of (15) is true if 
(k — A\APT™*)g + ((L— AAPH™1)(E +a) — la)g’ < 0. (18) 


By (16), (17), the inequality (18) is equivalent to 


—01 APT (E +a)” + abalE + a)?—* — O3(da)? < 0, 


with 
ğa 1 — Abtm} a (G+ 2)(q—- @-—m) 
B+m-1’ °? (@+m—1)((B+2)p—(28+m—))’ 
ee: aie B+1 
5 2BFm-1 (8+2)p—Q8+m-—1) 
Let 


(8 +1)(26 +m- 1) 


B+m—-1 
to A > ea - GB +m— 1)’ 
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then 6; > O(i = 1,2,3). Set p(y) = —@,A9+™-1y3=7 + abzy — 03(da)7 with y = (£ +a)7-}, 
then y(y) is a concave function, and attains its maximum at 
« _ {abaly -1)\ 7t 
s= 6 _ 
1y 
So, if we choose 
= >i 
d > aa yr 
O3y\ AY 
then y(y*) < 0, and thus (18) is true. 


Furthermore, if 
b2 baly = 1) JET 
Fa 
> max {1 go ( 01y ) i, 
with a small enough, we have 





(AC)P-8-™ (a7 — 1) 


=1)(B4+1) (6+2)p—(26+m+41 B+1 
epee) eee cee ene e ) , 


B+m-—-1 
which yields the second inequality of (15). 
In summary, @ is a supersolution of (1) with uo small enough if p > pe =28+m+1. 
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